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Research interest in the conversion of lignocellulosic biomass into energy and fuels through the pyrolysis 
process has increased significantly in the last decade as the necessity for a renewable source of carbon 
has become more evident. For optimal design of pyrolysis reactors, an understanding of the pyrolysis 
kinetics of lignocellulosic biomass is of fundamental importance. The distributed activation energy 
model (DAEM) has been usually used to describe the pyrolysis kinetics of lignocellulosic biomass. In this 
review, we start with the derivation of the DAEM. After an overview of the activation energy distribution 
and frequency factor in the DAEM, we focus on the numerical calculation and parameter estimation 
methods of the DAEM. Finally, this review summarizes recent results published in the literature for the 
application of the DAEM to the pyrolysis kinetics of lignocellulosic biomass. 


© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The pyrolysis of lignocellulosic biomass has attracted considerable 
attention because it is the thermochemical process that converts 


* Corresponding author at: Biomass Energy Engineering Center, Key Laboratory of 
Urban Agriculture (South) Ministry of Agriculture, School of Agriculture and 
Biology, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, 
PR China. Tel./fax: +86 21 34206624. 

E-mail address: jmcai@sjtu.edu.cn (J. Cai). 


http://dx.doi.org/10.1016/j.rser.2014.04.052 
1364-0321/© 2014 Elsevier Ltd. All rights reserved. 


lignocellulosic biomass into liquid fuel (bio-oil), which is easy to store 
and transport when compared to raw biomass [1,2]. Knowledge of the 
kinetics is important in the computational fluid dynamics modeling of 
biomass pyrolysis process [3]. And in the pyrolysis process of 
lignocellulosic biomass, the relative rates of decomposition, cracking, 
and repolymerization/condensation reactions influence the quantity 
and quality of bio-oil produced as well as the long-term storage 
stability of the bio-oil [4]. 

The distributed activation energy model (DAEM) is a multiple 
reaction model, which is widely used in the pyrolysis of 
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lignocellulosic biomass [5]. The model assumes the decomposition 
mechanism take a large number of independent, parallel, first- 
order or nth-order reactions with different activation energies 
reflecting variations in the bond strengths of species. The differ- 
ence in activation energies can be represented by a continuous 
distribution function [6]. The DAEM is not only used to describe 
the pyrolysis kinetics of lignocellulosic biomass and its main 
components, but also to describe the pyrolysis kinetics of coal 
[7-9] and other thermally degradable materials [10-13]. 

In order to better understand the DAEM and its application in 
the pyrolysis of lignocellulosic biomass, an overview of the 
research and development in the DAEM for the pyrolysis of 
lignocellulosic biomass is conducted. The derivation of the DAEM 
was first presented. Then, the activation energy distribution and 
frequency factor in the DAEM are overviewed and the numerical 
calculation of the DAEM is discussed. In Section 5, the methods for 
the estimation of the DAEM kinetic parameters are critically 
studied. Finally, the up-to-date research efforts in the DAEM for 
the pyrolysis of lignocellulosic biomass are summarized. 


2. Derivation of DAEM 


The concept of an activation energy distribution was originally 
proposed by Vand [14] and later used by Pitt [15] and Anthony [16]. 
The DAEM for lignocellulosic biomass pyrolysis may be applied to 
either the total amount of volatiles released, or to the amount of an 
individual volatile constituent. The description in this work follows 
the total amount volatiles released. It is assumed that lignocellulosic 
biomass has many chemical groups which are numbered i=1....,s, 
and the released mass fraction for the ith group is V(t). The 
contribution of the decomposition of the ith group is assumed to 
be described by the pseudo-nth-order rate equation 


dVi/Vi)_, (VtX-Vi\"_ a E;\ VŽ- VA" 
dt =k; Ve = Aj; exp -RT Ve a) 


where V;* is the total released mass fraction for the ith group, n is 
the reaction order and k; is the reaction rate constant, A; is the 
frequency factor, E; is the activation energy, R is the universal gas 
constant, T is the temperature, and t is the time. 

Eq. (1) can be written in its integral form 
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The complexity of lignocellulosic biomass pyrolysis is such that 
a continuous distribution f(E) of activation energies is assumed 
where JE ven f(E)dE describes the probability that chemical groups 
with a sample have an activation energy in this range E to E+ AE. 
The mass fraction of potential volatile material with activation 


energies between E and E+dE, dV*, at a given time t is V*f(E)dE. 


dV* = V*f (E)dE (3) 


Then, V;* and V; are replaced by dV* and dV, respectively. Then 
Eq. (2) becomes 
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The integration of Eq. (4) leads to 
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where a is the degree of conversion. 
Differentiating Eq. (5) with respect to t, the expression of 
da/dt is 
o0 n/(1—n) 
da | Jo A exp (—@) -a-m Joa exp (—)dt| f(EdE n#1 
dt Jo A exp |-4- fiA exp (—#)dt| fd n=1 
(6) 


Most laboratory experiments involving the pyrolysis of ligno- 
cellulosic biomass are performed at temperatures linearly varying 
with time from a starting temperature, To [17]. Under the linear 
heating program, Eqs. (5) and (6) can be rewritten as Eqs. (7) and (8) 
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where £ is the heating rate. 


3. fE) and A in DAEM 


A popular choice for the activation energy distribution is the 
Gaussian distribution centered at Eo with standard deviation o 
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The disadvantage of the Gaussian distribution is that it is 
symmetric whereas the actual reactivity distributions tend to be 
asymmetric [18]. The asymmetry can be accounted for by repla- 
cing the Gaussian distribution with other distributions, e.g. the 
Weibull [19-21] or Logistic [22,23] distributions. 

The Weibull distribution was first proposed to describe chemi- 
cal reactivity distributions for the pyrolysis of petroleum source 
rocks [19]. For the Weibull distribution, f(E) is given by 


A =e E 
Ey=—“(-—* — 10 
f(E) =( ; e| 7 (10) 


where 4 is the shape parameter, 7 is the width parameter, and y is 
the activation energy threshold (E>y). The mean activation 
energy Eo and the standard deviation o of the Weibull distribution 
are defined by 


1 
Eo=r+nr(1+5) (11) 
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where I is the Gamma function. Because of the activation energy 
threshold, Eqs. (7) and (8) should be written in the following forms 
if the Weibull distribution is considered 
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The Logistic distribution was presented in our previous paper 
[22] and was successfully applied to describe the pyrolysis kinetics 
of cellulose [23] and grape residues [24]. For the Logistic distribu- 
tion, f(E) is given by 


m exp [—(a(E-p)/V30)] 
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where „ and o are the mean value and standard deviation of the 
Logistic distribution. 

Capralis et al. studied the pyrolysis of two coals and found that 
the coal pyrolysis process can be divided into two steps: the tar 
and light hydrocarbon gas formation during the primary pyrolysis 
and the char condensation, cross-linking reactions, and a further 
gas production during the secondary pyrolysis [9]. To better 
describe the kinetic behavior, the authors proposed a double- 
Gaussian distribution to present the reactivity distribution [9] 


fB = wf (E)+ 1 —w)fa(E) (16) 


where f,(E) and f2(E) are Gaussian distributions, and w is the 
weight parameter. 

There are three ways to deal with the frequency factor in the 
literature. The conventional method to use the DAEM is to select a 
constant frequency factor value for all reactions (i.e. 1.67 x 10”? s~! 
for coal pyrolysis [8], and 10'413, 107! and 10'3°°s~' for the 
pyrolysis of hemicellulose, cellulose, and lignin, respectively [25]). 
If A is a constant, Eqs. (7) and (8) can be rewritten into the 
following equations 
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In the two above equations, the inner dT integral is 
T E 
i exp( =r) dT = KE,T) — KE, To) a9 
Ta RT 


where I(E,T)= R exp (—(E/RT))dT is the temperature integral 
[26,27]. 

Second, the frequency factor in the DAEM is assumed to be 
connected with the temperature 


A=AoI™ (20) 


where Ao is a constant and the exponent m ranges from — 1.5 to 
2.5. This dependence of the frequency factor on the temperature in 
the DAEM was found to be used to describe the pyrolysis kinetics 
of lignin [28], peanut shell, pistachito shell, and sunflower shell [29]. 


In this case, Eqs. (17) and (18) become 
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In Eqs. (21) and (22), the inner dT integral becomes 
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where Im(E, T) = i T™exp (—(E/RT))dT is the general temperature 
integral [30]. Here we especially mention the temperature integral 
or the general temperature integral in order to better state the 
numerical solution of the DAEM in the next section. 

Finally, the frequency factor in the DAEM is assumed to follow 
the compensation effect (the effect of an increase of Eo is partially 
or completely offset by a ‘compensatory’ increase in A [31]). 
The relationship between Eo and A is most usually represented 
[25,29,32] 


InA=a+bEo (24) 


where Eo is the mean value of activation energy distribution, a and 
b are constants. The value of a is usually close to zero [33]. 


4. Numerical calculation of DAEM 


Due to the fact that there is an inner dT integral and outer dE 
integral in the DAEM, it is very difficult to obtain the exact 
analytical solution of the DAEM. In the literature, many mathe- 
matical approximate approaches were pursued to deal with the 
DAEM. Suuberg expressed the temperature integral over E in terms 
of complementary error function [34]. Du et al. used the step 
function approximation technique to determine the activation 
energy distribution [35]. Donskoi and McElwain used high order 
Gauss—Hermite integration to evaluate the outer dE integral [36]. 
Rostami et al. presented how it is possible to obtain an approx- 
imate closed form of the model by assuming that as conversion 
proceeds, the functional groups with the lowest activation ener- 
gies are lost first [37]. 

Since it is difficult to analytically solve the DAEM, several 
researchers dealt with the DAEM by means of numerical techni- 
ques. Braun and Burnham replaced the upper and lower limits of 
the outer integral with Eo+30 and Eg—3o0 and divided the 
activation energy space into 96 equal intervals [38]. Güneş and 
Güneş used 500 kJ mol! as the upper limit of the outer dE 
integral and 50 integral intervals [39]. 

For purposes of numerical calculation of the DAEM, either 
general purposed mathematical software or a computer program 
developed in any programming language is used. In this work, we 
developed a program written in the Mathematica software system 
to numerically solve the DAEM. 

Either the temperature integral or the general temperature 
integral has no exact analytical solution [40], but they can be 
numerically solved. In the Mathematica software system, the 
temperature integral and general temperature integral can be 
expressed in the following forms 


I(E, T) = ; Gamma |- 1, Z (25) 
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E m+1 E 

Im(E, T) = (3) Gamma |- 1= maa (26) 

For the outer dE integral, we use the ‘NIntegrate’ function of the 
Mathematica software system with the integral limits 0 and 
Eo+ 300. 

To illustrate the proceeding method, the following theoretically 
simulated DAEM processes with Gaussian distribution are consid- 
ered. The kinetic parameter values of the first process are as 
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Fig. 1. The simulated a—T and da/dT—T curves of the Gaussian DAEM processes with 
the following parameters: A= 1.67 x 10”? s~ 1, Eọ=228 kj mol~!, =10 K min™ t, n=1.0, 
and g=0.5% Eo, 2.5% Eo, 5.0% Eo, 7.5% Eo and 10.0% Eo. 
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Fig. 2. The simulated a—T and da/dT—T curves of the Gaussian DAEM processes 
with the following parameters: A= 1.67 x 10° s~!, Ep =228 kJ mol~!, 6=10 K min™!, 
o=5.0% Eo, and n= 1.0, 1.5, 2.0, 2.5 and 3.0. 


Table 1 


follows: A=1.67 x 10” s~', Eo=228 kJ mol~', 6=10 K min~', n=1 
and ¢=0.5%Eo, 2.5% Eo, 5.0% Eo, 7.5% Eg and 10.0% Eo. The values of 
A, Eo and £ were taken from the literature [8], which were used to 
describe the pyrolysis kinetics of coal. The parameter values in the 
second process are identical to the first example except for 
o0=5.0% Eo and n=1.0, 1.5, 2.0, 2.5, and 3.0 (those n values were 
used to describe the pyrolysis kinetics of biomass in the literature 
[41]). Figs. 1 and 2 show the a—T and da/dT—T curves for the first 
and second processes. From Figs. 1 and 2, it can be seen that 
increasing the values of o makes the reaction profile broader and 
more symmetric and that increasing the reaction order makes 
the reaction profile broader and makes it more skewed to high 
temperature. 


5. Parameter estimation methods 


The DAEM has some mathematical difficulties in parameter 
estimation because of the complex structure of the DAEM equation 
[28]. The parameter estimation methods for the DAEM can be 
classified as distribution-free and distribution-fitting methods. 
Miura and his colleague [42,43] proposed the differential and 
integral methods to estimate fE) in the DAEM without previous 
assumptions for E) in 1995 and 1998, respectively. In further 
discussion, they are referred to as the Miura differential method 
and the Miura—Maki integral method [6], which are distribution- 
free methods. The previous assumption for f(E) is required in the 
use of the distribution-fitting methods, which are based on some 
direct search optimization methods. 


5.1. Distribution-free methods 


The assumptions and detailed derivation processes of the 
Miura differential method and the Miura—Maki integral method 
can be found in the literature [42,43]. The procedures for calculat- 
ing f(E) by using these distribution-free methods are compared in 
Table 1. From Table 1, it can be obtained that the a—T data at 
different heating rates are required in the use of the Miura—Maki 
integral method, while the Miura differential method employs 
instantaneous a—T and da/dT—T data. Therefore, for the Miura 
differential method, it is more sensitive to experimental noise. 

To evaluate the validity of the Miura differential method and 
the Miura-Maki integral method, it is examine to see if fE) 
determined from the a—T and da/dT—T data constructed from 
known sets of DAEM parameters is consistent with the assumed 


Table 2 
DAEM parameters of cellulose, hemicellulose and lignin presented in the literature [25]. 


Material Eo (kj mol~!) o (kJ mol~!) A 
Hemicellulose 175.6 4.3 101452 
Cellulose 185.4 3.5 102° 
Lignin 195.4 32.0 1013-98 


Comparison between the Miura differential method and the Miura—Maki integral method to estimate f(E). 


Step The Miura differential method 


The Miura-Maki integral method 


1 Obtain a versus T and da/dT versus T relationships at three heating rates at least 


Calculate the values of In[(da/dT)/(1 —a)] at selected a values for different heating 


rates 
Plot In[{(da/dT)/(1 —a)] versus 1/T at selected a values, and determine E 


Ww 


Obtain a versus T relationships at three heating rates at least 

Calculate the values of In(,/T”) at selected a values for different heating 
rates 

Plot In(4/T?) versus 1/T at selecteda values, and determine E 

Plot a versus E, and differentiate a versus E relationship to obtain f(E) 


4 Plot a versus E, and differentiate a versus E relationship to obtain fE) 
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Fig. 3. The activation energy values as a function of conversion for the pyrolysis of hemicellulose, cellulose and lignin. 
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Table 3 
Different sets of kinetic parameters of the first-order Gaussian DAEM for cellulose 
pyrolysis. 


Set of kinetic Alst Eo o SSR 
parameters (kJmol~!) (kJ mol~') 
1 5.5472 x 10'! 161.8438 4.3217 0.012030 
2 6.3165 x 10" 185.1297 4.5979 0.002649 
3 1.0107 x 10'* 187.4462 4.5818 0.002569 
4 9.9504 x 10'4 198.7059 5.2204 0.003468 
5 1.0004 x 10'® 210.0828 5.8349 0.004290 
6 9.8514 x 107 232.7228 6.9938 0.005719 
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one. For this purpose, some simulated DAEM processes are 
considered. 

In the literature [25], the first-order Gaussian DAEM was used 
to successfully describe the pyrolysis kinetics of hemicellulose, 
cellulose and lignin, which are three main components of ligno- 
cellulosic biomass. Their parameters are employed to construct the 
simulated DAEM processes, as shown in Table 2. 

Fig. 3 shows the activation energy values as a function of 
conversion for the pyrolysis of hemicellulose, cellulose and lignin. 
From Fig. 3, it can be observed that the activation energies 
significantly increase with increasing conversion. Differentiating 
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Fig. 5. Multiplicity in kinetic parameters of the first-order Gaussian DAEM for cellulose pyrolysis. 
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the a versus E relationships gives f(E) curves, which are shown in 
Fig. 4 with the pre-assumed Gaussian distributions. The estimated 
activation energy distributions are close with the assumed dis- 
tribution for lignin. However, the errors involved in the activation 
energy distributions obtained from the Miura differential method 
and the Miura—Maki integral method are significant for hemicel- 
lulose and cellulose. Therefore, the validity of the Miura differ- 
ential method and the Miura-Maki integral method for the 
estimation of the DAEM kinetic parameters should be system- 
atically investigated. 


5.2. Distribution-fitting methods 


The conventional procedures of the distribution-fitting meth- 
ods for the determination of the kinetic parameters of DAEM 
include four steps 


(1) A is assumed to be an unknown constant, or is directly 
assigned to a given constant, or A is equaled to AoT”; 
(2) KE) is assumed to be a Gaussian or another distribution; 
(3) The following objective function is constructed 
Na 
SSR = 2 (Qexp,i — acai) (27) 
t= 
where aexp and acai are experimental and calculated values of 
conversion, respectively; 
(4) A and fE) are obtained by minimizing the objective function (27). 


In the literature, the grid search method [19,44], genetic algo- 
rithm [45], pattern search method [46], Multistart algorithm [47], 
Hook-Jeeves method [48], and simulated annealing optimization 
method [28] were usually used for minimizing Eq. (27). 

The disadvantage of the distribution-fitting methods is that 
different sets of DAEM kinetic parameters provide an approxi- 
mately equally good fit to experimental data at a single heating 
rate [25,49]. To demonstrate the existence of multiplicity in DAEM 
kinetic parameters, the experimental data (52 data points) of 
cellulose pyrolysis at a heating rate of 10K min~! is used. The 
Gaussian distribution of first-order reactions is employed to 
analyze the experimental data of cellulose pyrolysis. Parameter 
estimation is performed by means of the pattern search method. 
Several sets of parameters is obtained with different initial guesses 
and listed in Table 3. The results from the first-order Gaussian 
DAEMs with different sets of kinetic parameters are shown in 
Fig. 5, where it can be seen that the DAEMs compare very well 
with the experimental data. Fig. 6 shows the relationship between 
the multiple sets of the mean value (Eg) of the activation energy 
distribution and the corresponding frequency factor (A). A perfect 
linear relationship can be obtained between InA and Eo 


In A = 28.61335 + 4.92587E0 (28) 


where A and Ep are expressed in s~' and kJ mol~', respectively. 


This kind multiplicity of the DAEM kinetic parameters was studied 
in detail in the literature [12,49]. Therefore, the use of the method 
involving nonlinear least squares optimization of multiple heating 
rate data is expected to result in unique kinetic parameters [50]. 


6. DAEM for lignocellulosic biomass pyrolysis kinetics 


Lignocellulosic biomass consists of three major components 
(i.e. cellulose, hemicellulose, and lignin) and also contains a variety 
of minor components (i.e. inorganic matters) [51]. Cellulose, 
hemicellulose, and lignin are all long-chain biopolymers [52,53]. 
There are a lot of reactions involved in the thermal decomposition 
of cellulose, hemicellulose, and lignin. Therefore, it is not accurate 
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Fig. 6. The relationship between Eo and InA for multiple sets of kinetic parameters. 
The solid line represents the linear fitting. (A is expressed in s~'). 


enough to describe the entire pyrolysis processes of them with a 
single reaction [54]. In fact, the pyrolysis of a biomass component 
requires a range of activation energies that vary with the size of the 
polymer and the location of the bonds that are being broken [55]. 
Therefore, the DAEM has been used previously to describe the 
pyrolysis of biomass components [23,55-60]. 

In general, the conversion rate curve (i.e. the da/dT— T curve) of the 
pyrolysis of lignocellulosic biomass shows a peak and a peak shoulder 
on the left side of the peak [5]. However, according to the numerical 
results of the DAEM processes [29,39,41,61,62], the conversion rate 
curves of the DAEM with a certain activation energy distribution (i.e. 
Gaussian, Weibull, Logistic distributions) and certain frequency factor 
(i.e. A=constant, A=Ag7™, A = e*+ £0) show one single peak and no 
peak shoulder. Thus, it is impossible to describe the entire pyrolysis 
process of lignocellulosic biomass with such the DAEM. Some 
researchers reported that they used the DAEM to describe the 
pyrolysis kinetics of lignocellulosic biomass based on the Miura—Maki 
integral method [63-67]. However, they only obtained some values of 
E, and A, or the ranges of E and A in a certain range of conversion. fE) 
was supposed to be given by numerical differentiating the a versus Ex 
relationship. However, sometimes there were several points in their 
obtained a—E, plot, sometimes a varied with E, dramatically. In those 
situations, it is difficult to get the derivatives of the a versus E, 
relationship. 

Since the pyrolysis kinetics of biomass components can be 
described well by using the single DAEM, some comprehensive 
models consisting of two or three parallel DAEM reactions were 
proposed to describe the overall decomposition kinetics of lignocellu- 
losic biomass [25,54,55,68]. In those models, the lignocellulosic 
biomass was regarded as a sum of several pseudocomponents. It 
was further assumed that there was no interaction among these 
pseudocomponents and the pyrolysis of each pseudocomponent 
would be described by using a single first-order Gaussian DAEM. 


M PeR 
a(t) = di =e) (29) 
i=1đi 
da E 1Ci(da;/dt) 
da yw, (30) 
a(t) =1-—— [ex - [mex (i)i EE dE 
O 2a Jo p a TSENTRI 20;2 
(31) 


Table 4 
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Summary of DAEM for lignocellulosic biomass pyrolysis kinetics. 


Biomass 


Jerusalem 
artichoke tubers 


Pine pellets 


Sugarcane trash 
sugarcane 
bagasse pine 
sawdust 

Cedar sawdust 
coffee residue 


Xylan 
Microcrystalline 
cellulose 


Indulin AT lignin 
Lignoboost™ 
lignin 


Acetocell lignin 


Alcell lignin 


Kraft lignin 


Alcell lignin 


Hemicellulose 
Cellulose 
Lignin 

Wood 


Fresh water algae 


Experimental method and condition 


Thermogravimetric analysis—linear heating: 


5, 10, 20 and 30 K min~!; inert atmosphere: 


nitrogen 

Thermogravimetric analysis—linear heating: 
10, 15, and 20 K min!; inert atmosphere: 
nitrogen 


Model 


The first-order DAEM 


The first-order DAEM 


Thermogravimetric analysis (sugarcane trash, The first-order DAEM 


sugarcane bagasse, pine sawdust, cedar 
sawdust, coffee residue)—linear heating: 
0.0833-1 K s~!; inert atmosphere: nitrogen 


Thermogravimetric analysis—linear heating: 
5 K min~!; inert atmosphere: helium 


Thermogravimetric analysis—linear heating: 


5,10, and 15 K min~! for indulin AT lignin, 5, 


15, and 40 K min”! for Lignoboost™ and 
acetocell lignin; inert atmosphere: nitrogen 


Thermogravimetric analysis—linear heating: 
5,10 and 15 K min” !; inert atmosphere: 
nitrogen 


Fixed-bed reactor—linear heating rates: 5, 10 
and 15 K min” !; inert atmosphere: helium 


Thermogravimetric analysis—linear heating: 
5, 10, and 15 K min~!; inert atmosphere: 
nitrogen 


Thermogravimetric analysis—linear heating: 


5 K min” !; inert atmosphere: nitrogen 


Thermogravimetric analysis—linear heating: 
5 and 10 K min~!; inter atmosphere 


The first-order Gaussian 


DAEM 


The first-order Gaussian DAEM 


The first-order DAEM 


The first-order Gaussian DAEM with 


A=AgI™ 


The first-order Gaussian DAEM 


The nth-order Gaussian DAEM 


Parameter 
estimation method 


The Miura—Maki 
integral method 


The Miura-—Maki 
integral method 


The Miura-Maki 
integral method 


Distribution-fitting 
method: pattern 
search method 


The Miura-Maki 
integral method 


Miura—Maki 
integral method 


Distribution-fitting 
method: simulated 
annealing 
optimization 
method 
Distribution-fitting 
method 


Distribution-fitting 
method: Multistart 
algorithm 


Kinetic parameters 


E: 152-232.5 kJ mol! 


A: 2.93 x 101 

7.97 x 10” s7! 

E: 160-270 kj mol~? 
A: 10''-10' s~! 


Sugarcane trash 
E: 110-320 kJ mol`! 
Sugarcane bagasse 


E: 110-320 kJ mol~! 
Pine sawdust 

E: 200-350 kJ mol~! 
Cedar sawdust 

E: 120-315 kJ mol~! 
Xylan— 

E9=178.311 kj mol` t, 
o=5.848 kJ mol~!, 
A=10!2:702 s7 T 
Cellulose— 
Eo=210.037 kJ mol~, 
o=0.944 kJ mol` t, 
A= 1013-904 s7 q: 
Indulin AT lignin: 
Eo=192 kJ mol~!, 
o=43 kJ mol~!, 
A=10% s7! 
Lignoboost™ lignin: 
Eo=193 kJ mol~!, 
o=37 kj mol`’, 
A=10” s7! 
Acetocell ligni 
Eo=193 kJ mol~!, 
o=47 kj mol` t, 
A=10” s7! 

Alcell lignin—E: 129- 
361 kJ mol~', A: 


= 


6.2 x 10"'-9,2 x 107? s71, 


Fo=254.3 kJ mol!, 
o=29.0 kJ mol~! 
Kraft lignin—E: 80- 
158 kj mol” t, A: 


3.3 x 10’-1.84 x 10° s71, 


E=123.0 kJ mol~!, 
o=11.0 kJ mol`! 
Eo=135 kJ mol~! 
o=33kj mol"! 
m=2.8 

Ao=92670 min~! K~! 


Hemicellulose— 
Eo=132.9 kJ mol}, 
o=24kj mol™ t, 
A=1.25 x 101° s~! 
Cellulose— 
Eo=175.6 kJ mol~!, 
o=04kj mol~ t, 
A=3.41 x 10”? s7! 
Lignin— 

Eo=101.0 kJ mol}, 
o=11.3 kJ mol` t, 
A=2.22 x 10°s~! 
Wood— 

E=123.52 kj molt, 
o=5.97 kJ mol” !, 
A=1.57 x 10° s7! 
5K min™!: 
Eo=189.15 kJ mol~!, 
o=14.73 kJ mol` t, 
A=11291 x 10" s~?, 
n=7.88 

10 K min™!: 
Eo=190.02 kj mol ~’, 
o=14.73 kj mol` t, 
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References 
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[55] 


[58] 


[57] 


[28] 


[56] 
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Table 4 (continued ) 


Biomass Experimental method and condition Model Parameter Kinetic parameters References 
estimation method 
A=2.6021 x 10'° s71, 
n=6.63 
Cellulose Thermogravimetric analysis—linear heating: The nth-order Logistic DAEM Distribution-fitting =258.5718 kJ mol`! [23] 
5, 25, and 50 K min~!; inert atmosphere: method: pattern o=2.6601 kJ mol7! 
nitrogen search method A=1.6218 x 10” s~! 
n=1.1101 
Corn stalk skin Thermogravimetric analysis—Linear heating: The first-order DAEM The Miura-Maki E: 62-169 kJ mol`! [65] 
20, 35, 50, 75, and 100 K min™ t; Inert integral method A: 10108-1065 s~! 
atmosphere: nitrogen 
Rice straw Thermogravimetric analysis—linear heating: The first-order DAEM The Miura-Maki Rice straw—E: 118- [64] 
Rice husk 2, 5 and 10 K min” !; inert atmosphere: integral method 208 kJ mol~! 
Corncob nitrogen fE) peaks at 
170 kJ mol`! 
Rice husk—E: 150- 
224 kJ mol~! 
fE) peaks at 
174 kJ mol`! 
Corncob—E: 140- 
229 kJ mol`! 
fE) peaks at 
183 kJ mol~! 
Sawdust Thermogravimetric analysis—linear heating: The first-order DAEM The Miura-Maki E: 161.9-202.3 kJ mol~! [69] 
5, 10, 15, and 20 K min~!; atmosphere: integral method A: 4.8 x 10°- 
hydrogen or syngas (H3: 68.1%, CO: 30.1%, C4- 4.5 x 10”? s7! 
C4: 0.9%, CO2: 0.9%) 
Wheat straw Thermogravimetric analysis—linear heating: Two parallel first-order Gaussian Distribution-fitting See Table 4.1 [68] 
Winter barley 11, 22, and 47 K min` !; inert atmosphere: DAEM reactions method 
straw nitrogen 
Oat straw —dm/dt = cı daı /dt+c2daz /dt 
where a and az are described by 
the first-order Gaussian DAEM 
Corn stalk Thermogravimetric analysis—linear heating: Three parallel first-order DAEM Distribution-fitting See Table 4.2 [25] 
Rice husk 4 and 40 K min“! (corn stalk, rice husk, reactions method 
Sorghum straw sorghum straw, and wheat straw); stepwise 
Wheat straw program: 20 K min! with steps at 250, 300, —dm/dt = 3_ ,cjda;/dt 
350, and 400 °C, 20 K min~! with steps at 225, where c; is the amount of volatiles 
275, 325, and 375 °C for corn stalk, sorghum formed from a unit mass of 
straw and wheat straw; 20 K min~! with pseudocomponent j; aj is described 
steps at 225, 275, 325, 375, and 425 °C, by the first-order Gaussian DAEM. 
20 K min~! with steps at 200, 250, 300, 350, 
and 400 °C for rice husk; Inert atmosphere: 
argon 
Corn stover Thermogravimetric analysis—linear heating: Three parallel first-order DAEM Distribution-fitting See Table 4.3 [55] 


Cotton stalk 

Palm oil husk 
Pine wood 

Red oak 
Sugarcane bagasse 
Switchgrass 
Wheat straw 


5 and 20 K min` !; inert atmosphere: helium reactions 


method: pattern 
search method 


—dm/dt = ¥3_ ,cjda;/dt 


where cj is the amount of volatiles 
formed from a unit mass of 
pseudocomponent j; a; is described 


by the first-order Gaussian DAEM 


Table 4.1 
Kinetic parameters of lignocellulosic biomass pyrolysis in the literature [68]. 


Wheat straw Winter barley straw Oat straw 
Eo: (kJ mol~') 167.3 167.3 167.3 
Eo2 (kJ mol~') 225.7 225.7 225.7 
c (kJ mol!) 1.8 1.8 1.8 
02 (kJ mol~') 25.9 25.9 25.9 
Ay (s7 t 101253 101291 101221 
A2 (s~ 1) 1017-78 101814 1017-20 
Cy 0.36 0.31 0.26 
C2 0.35 0.36 0.47 


where the weight factor c; is equal to the amount of volatiles formed 
from the ith pseudocomponent, and M is the number of 
pseudocomponents. 

It is noted that a pseudocomponent represents the totality of 
those decomposing species that can be described by the same 


kinetic parameters. Those pseudocomponents are usually linked 
the real biomass components. In particular, for the three-parallel- 
DAEM-reaction model [25,55], the first and second pseudocompo- 
nents represent the fractions of hemicellulose and cellulose that 
react at low temperature, and the third pseudocomponent repre- 
sents the sum of the remaining amounts of hemicellulose and 
cellulose plus the fractions of lignin. 

By following the above stepwise discussion, the research and 
development efforts in the literature on the DAEM for lignocellu- 
losic biomass pyrolysis are summarized in Table 4. Most kinetic 
studies of lignocellulosic biomass pyrolysis have been made using 
thermogravimetric analysis (TGA) at relatively slow heating rates. 
In the tabulated summary, the elemental composition of lignocel- 
lulosic biomass sample, the experimental methods and conditions 
are highlighted. The resultant DAEM kinetic parameters together 
with the parameter estimation methods are also provided. The 
efforts on DAEM for the pyrolysis of lignocellulosic biomass 
indicate that the model consisting of three parallel DAEM 


Table 4.2 
Kinetic parameters of lignocellulosic biomass pyrolysis in the literature [25]. 
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Corn stalk Sorghum straw Rice husk Wheat straw 

Eo, (kJ mol~!) 176 176 176 176 

Eo2 (kJ mol~') 185 185 185 185 

Eos (kJ mol ') 189 189 189 189 

cı (kJ mol`!) 71 71 71 71 

G (kJ mol~') 17 {7 17 1.7 

3 (kJ mol~') 32.7 32.7 32.7 32.7 

Ay (s~!) 10143 10143 10143 10143 

Az (s~!) 1013-71 101371 1071 101371 

A3 (s~!) 1013-99 1013-99 1013-99 1013-99 

C1 0.14 0.10 0.08 0.16 

C2 0.25 0.23 0.36 0.38 

C3 0.31 0.34 0.33 0.19 

Table 4.3 
Kinetic parameters of lignocellulosic biomass pyrolysis in the literature [55]. 
Corn stover Cotton stalk Palm oil husk Pine wood Red oak Sugarcane bagasse Switchgrass Wheat straw 

Eo, (kJ mol!) 179.602 178.188 169.710 186.701 183.109 184.750 186.776 175.506 
Eo2 (kJ mol~') 207.381 205.287 199.966 214.364 209.58 212.483 212.057 204.244 
Eos (kJ mol!) 239.344 239.463 236.110 271.758 242.145 234.759 260.952 240.610 
cı (kj mol!) 5.888 5.591 5.147 8.769 8.342 5.419 4.586 5.375 
Ga (kJ mol~') 1.174 1.789 1.313 1.126 0.706 1.339 1.361 1.130 
63 (kJ mol™') 31.414 41.767 39.998 29.250 26.583 36.493 39.293 38.422 
Ay (s~!) 10 13.007 1013-132 1012-999 1013-006 1013-207 1013-007 1013-007 1012700 
A2 BH 1013-949 1013-800 1013618 1013-803 1013-618 1913-805 1013-804 1013-800 
A3 E 1016-003 1015-738 1015-953 1016399 10 15.006 1015-769 1016536 1015812 
c 0.16337 0.12999 0.10116 0.32147 0.27282 0.22009 0.23651 0.15916 
C2 0.25386 0.35824 0.32121 0.38557 0.37320 0.35035 0.33376 0.28027 
C3 0.16481 0.20008 0.19241 0.08412 0.12204 0.12289 0.20159 0.23410 


reactions is the most accurate and up-to-date approach for 
modeling the pyrolysis kinetics of lignocellulosic biomass. 


7. Summaries 


The pyrolysis of lignocellulosic biomass involves many reac- 
tions, which can't be well described by a simple single kinetic 
model. In the literature, the DAEM was usually used to describe its 
pyrolysis kinetics. This review presents a state-of-the-art review of 
the DAEM for lignocellulosic biomass pyrolysis. This review starts 
with the derivation of the DAEM and a brief overview of the 
activation energy distribution and frequency factor in the DAEM. 
Then the parameter estimation methods for the determination of 
the DAEM kinetic parameters are critically studied. The final part is 
the overview of the application of the DAEM to lignocellulosic 
biomass pyrolysis. A handy, up-to-date overview of the existing 
research and developments efforts and major findings in the 
DAEM for lignocellulosic biomass pyrolysis is provided. It is hoped 
that this review is critical for helping understand the application 
of the DAEM in the thermal decomposition of other degradable 
materials. 
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